Properties of nitrogen-vacancy centers in diamond: group theoretic approach 
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We present a procedure that makes use of group theory to analyze and predict the main 
properties of the negatively charged nitrogen-vacancy (NV) center in diamond. We focus 
on the relatively low temperatures limit where both the spin-spin and spin-orbit effects are 
important to consider. We demonstrate that group theory may be used to clarify several 
aspects of the NV structure, such as ordering of the singlets in the (e 2 ) electronic configu- 
ration, the spin-spin and the spin-orbit interactions in the (ae) electronic configuration. We 
also discuss how the optical selection rules and the response of the center to electric field 
can be used for spin-photon entanglement schemes. Our general formalism is applicable to 
a broad class of local defects in solids. The present results have important implications for 
applications in quantum information science and nanomagnetometry. 



I. INTRODUCTION 



During the past few years nitrogen- vacancy (NV) centers have emerged as promising candidates 
for a number of applications [IH1] ranging from high spatial resolution imaging [5j to quantum 
computation [6j. At low temperatures, the optical transitions of the NV center become very 
narrow and can be coherently manipulated, allowing for spin-photon entanglement generation [7] 
for quantum communication and all optical control [8]. A detailed understanding of the properties 
of this defect is critical for many of these applications. Several studies have addressed this issue 
both experimentally [HI [10] and theoretically I12j. Furthermore, other atom-like defects can 
potentially be engineered in diamond [13] and other materials with similar or perhaps better 
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properties suitable for the desired application. Therefore, it is of immediate importance to develop 
a formalism to analyze and predict the main properties of defects in solids. 

Here we present a formalism based on a group theoretical description. While we focus on 
describing the nitrogen-vacancy center in diamond, our formalism can be applied to any point 
defect in solid state physics. Our method takes advantage of the symmetry of the states to properly 
treat the relevant interactions and their symmetries. We apply group theory to find out not only 
the symmetry of the eigenstates but also their explicit form in terms of orbital and spin degrees of 
freedom. We show that this is essential to build an accurate model of the NV center. In particular, 
we analyze the effect of the Coulomb interaction and predict that the ordering of the triplet and 
singlet states in the ground state configuration is { 3 A2, 1 E, 1 Ai} and that the distance between 
them is on the order of the exchange term of the electron-electron Coulomb energy. This ordering 
has been debated over the last few years and our results agree with recent ah initio calculations 
carried in bulk diamond [14J. 

Our method is also used to analyze important properties of the center such as polarization 
selection rules. The explicit form of the states allows us to identify a particularly useful lambda- 
type transition that was recently used for spin-photon entanglement generation [7] • We also consider 
perturbations that lower the symmetry of a point defect, such as strain and electric field and how 
they affect the polarization properties. We also show that the non-axial spin-orbit interaction 
discussed in Ref. [15] does not mix the eigenstates of the center in a given multiplet. Instead, we 
find that the electron spin-spin interaction is responsible for the spin state mixing of the excited 
state as a result of the lack of inversion symmetry of the center. Finally, we analyze the effect 
of electric fields via the inverse piezoelectric effect and compare our results with experimental 
observations. We show that this effect can be used to tune the polarization properties of optical 
transitions and the wavelength of emitted photons, which is of direct importance for photon-based 
quantum communication between NV centers. Our study clarifies important properties of NV 
centers and provides the foundation for coherent interaction between electronic spins and photons 
in solid state. 

Our manuscript is organized as follows. In Section [TT] we present a general group theoretical 
formalism to calculate the electronic or hole representation of a point defect for a given crystal 
field symmetry and number of electrons contained in the defect. Next, we use group theory and 
the explicit form of the states to analyze the effect of the Coulomb interaction between electrons 
(Section |III[ ) and spin-spin and spin-orbit interactions for the NV center (Sections [V] and IV 



respectively). Next, we analyze the selection rules of the unperturbed defect in Section [Vl} Finally, 



in Section VII we analyze the effect of strain and electric field perturbations. 



II. STATE REPRESENTATION 



We are particularly interested in quasi-static properties of defects in crystals where the complex 
electronic structure can be observed spectroscopically. In this limit one can apply the Born- 
Oppenheimer approximation to separate the many-body system of electrons and nuclei. This 
approximation relies on the fact that nuclei are much slower than electrons. In this approximation 
the nuclei are represented by their coordinates and the physical quantities of the electrons depend on 
these coordinates as (external) fixed parameters. A defect in a crystal breaks down the translational 
symmetry reducing the symmetry of the crystal to rotations and reflections. These symmetries 
form a point group which in general is a subgroup of the point group of the lattice. The loss 
of translational symmetry indicates that the Bloch-states are no longer a good approximation to 
describe the point defect. In fact, some states can be very well localized near the point defect. 
These defect states are particularly important in semiconductors and insulators when they appear 
within the fundamental band gap of the crystal. 

In the tight binding picture, the electron system of the diamond crystal may be described as 
the sum of covalent-type interactions between the valence electrons of two nearest neighbor atoms. 
When defects involve vacancies, the absence of an ion will break bonds in the crystal, producing 
unpaired electrons or dangling bonds, cr,, which to leading order can be used to represent the single 
electron orbitals around the defect. The particular combination of dangling bonds that form the 
single electron orbitals {tp r } is set by the crystal field of the defect and can be readily calculated 
by projecting the dangling bonds on each irreducible representation (IR) of the point group of the 
defect [16J, 

Vr = P {r) cr l = l j i J2x i I ) Rea l , (1) 

where is the projective operator to the IR r, Xe is the character of operation R e (element) 
for the IR r, l r is the dimension of the IR r, and h is the order of the group (number of elements). 
A detailed application of Eq. ([I]) for the case of the NV center can be found in Appendix |Aj There 
are two non-degenerate totally symmetric orbitals «i(l) ad 02(1) that transform according to the 
one-dimensional IR Ai, and there are two degenerate states {e x ,e y } that transform according to 
the two-dimensional IR E. At this stage, group theory does not predict the energy order of these 
states. However a simple model of the electron-ion Coulomb interaction can be used to qualitatively 



4 



obtain the ordering of the levels [T7j . In Appendix [A] we model the effect of this interaction on 
the single electron orbitals, ip r , for the case of the NV center and find that the ordering of the 
states (increasing in energy) is ai(l), a\(2) and {e x , e y }. Indeed, ab initio density functional theory 
(DFT) calculations revealed [18, 19J that the a\(l) and a\{2) levels fall lower than the e x and e y 
levels, which demonstrates the strength of group theory for qualitative predictions. 

Once the symmetry and degeneracy of the orbitals are determined, the dynamics of the defect 
is set by the number of electrons available to occupy the orbitals. The orbitals with higher energy 
will predominantly set the properties of the defect. The spin character of the defect will be 
determined by the degeneracy of the orbitals and the number of electrons in them, leading to net 
spins S = {0, 1, 2, ...} if this number is even and S = ||, |, ... j if odd. 

In the case of the negatively charged NV center, each carbon atom contributes one electron, the 
nitrogen (as a donor in diamond) contributes two electrons, and an extra electron comes from the 
environment [19J, possibly given by substitutional nitrogens [20J. The ground state configuration 
consists of four electrons occupying the totally symmetric states and the remaining two electrons 
pairing up in the {e x ,e y } orbitals. In this single particle picture, the excited state configuration 
can be approximated as one electron being promoted from the ai(2) orbital to the e xy orbitals 

If two more electrons were added to any of these configurations, the wavefunction of the defect 
would be a singlet with a totally symmetric spatial wavefunction, equivalent to the state of an 
atom with a filled shell |21|, I22j . Therefore, the electronic configuration of this defect can be 
modeled by two holes occupying the orbitals e x , y in the ground state (e 2 electronic configuration) 
and one hole each in the orbitals a±{2) and e xy for the excited state (ae electronic configuration). 
A third electronic configuration, a 2 , can be envisioned by promoting the remaining electron from 
the orbital a±(2) to the orbitals e x ^ y . Hole and electron representations are totally equivalent and 
it is convenient to choose the representation containing the smallest number of particles. If a hole 
representation is chosen, some care must be taken, as some interactions reverse their sign, such 
as the spin-orbit interaction |21j. In what follows, we choose a hole representation containing two 
particles (instead of an electron representation containing four particles), since it is more convenient 
to describe the physics of the NV center. However, the analysis can be applied to electrons as well. 

The representation of the total n-electron wavefunction, including space and spin degrees of 
freedom, is given by the direct product of the representation of each hole Th n and its spin Tq, = 
Il n (j*hn ® Dij, where Di is the representation for a spin \ particle in the corresponding point 
group. The reduction or block diagonalization of the representation r<j, gives the eigenstates of the 
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hamiltonian associated with the crystal field potential and any interaction that remains invariant 
under the elements of the point group in question. These interactions include spin-orbit, spin-spin 
and Coulomb interactions, as well as expansions, contractions, and stress where their axes coincide 
with the symmetry axis of the defect. The eigenstates can be found by projecting any combination 
of the two electron wavefunction onto the irreducible representations of the group [16} [23] , 

^ = PrtipiW = \Y. xP*ReVlRe<P2, (2) 
e 

where ipi can be any of the orbitals in Eq. ([I]) and the subindex i refers to the hole i. In the 
case of the NV center, it is illustrative to note that the spin representation for the two particles 
can be reduced to D 1 / 2 <£> D 1 / 2 = A\ + A 2 + E, where A\ corresponds to the singlet state, and 
A 2 and E to the triplet state with zero and non-zero spin projections, respectively. A list of the 
eigenstates and their symmetries for the two hole representation can be found in Table |T] for the 
ground state (e 2 ) and the excited state (ae). For completeness, we include the doubly excited state 
(a 2 ) electronic configuration although this state is not optically accessible in the excitation process 
of the NV center in experiments. Note that each electronic configuration might have singlet and 
triplet states. The calculation performed to obtain Table [I] is similar to the calculation made to 
find the eigenstates when two spin particles are considered. However, in this case one should use 
the Wigner coefficients of the corresponding irreducible representation of the point group under 
consideration. 

Group theory can predict why the hyperfine interaction with the nuclear spin of the nitrogen 
in the excited state is more than an order of magnitude larger than in the ground state for both 
nitrogen species: the non-zero spin density in the ground state wavefunction of the NV center 
is mostly concentrated in the orbitals e xy , which have no overlap with the nitrogen atom. On 
the other hand, in the excited state, when one electron is promoted from the ai(2) orbital to one 
of the e x ,y orbitals, the non-zero spin density comes now from unpaired electrons occupying the 
orbitals <zi(2) and e x ,y As the orbital a\(2) is partially localized on the nitrogen atom, a sizable 
contact term interaction between the electronic spin and the nuclear spin of the nitrogen is expected 

PS El [22]. 

Up to now, eigenstates inside a given electronic configuration have the same energy, but the 
inclusion of the electron-electron Coulomb interaction will lift the degeneracy between triplets and 
singlets. The resulting energy splitting can be of the order of a fraction of an eV and it is analyzed 



for the ground state configuration of the NV center in Section III Furthermore, the degeneracy 
of triplet states is lifted by spin-orbit and spin-spin interactions of the order of GHz, where the 
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crystal field plays an important role. These interactions will be treated in sections IV andfV 



TABLE I: Partner functions of each IR for the direct product of two holes. The first column shows the 
electronic configuration and in parenthesis their triplet (T) or singlet (S) character. The last column shows 
the name of the state given in this paper and their symmetry. a(/3) stands for f (!) and E± — \ae± — e±a), 
where e± = T(e x ± ie y ), \X) = (|£_) - \E+))/2 and \Y) = (|£_) + \E+))i/2. 
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III. ORDERING OF SINGLET STATES 

For a given electronic configuration, the most relevant interaction is the electron-electron 
Coulomb interaction, which is minimized when electrons are configured in an antisymmetric spatial 
configuration. As the total wavefunction must be antisymmetric for fermionic particles, the spin 
configuration must be symmetric. As a result, the state with the largest multiplicity lies lower in 
energy. This analysis, known as the first Hund's rule, predicts that the ground state of the NV 
center should be the triplet S A2 state. We now address the question related to the order of singlets 
in the ground state electronic configuration e 2 . The order of singlet states has a great significance 
in understanding the spin-flipping fluoresence of the NV center, and ab initio DFT calculations 
were unable to address this issue properly due to the many-body singlet states. Since we have the 



7 



explicit form of the wavefunctions, we can work out the ordering of the singlets in a given electronic 
configuration by analyzing the expectation value of the Coulomb interaction, which can be written 
in the general form, 

C ab cd = J dV 1 dV2a*{r 1 )b*{r 2 )V(\r 1 -r 2 \)c(r l )d{r 2 ). 

Using this expression we find that in the ground state electronic configuration (e 2 ), the Coulomb 
interactions for these states are 

C{ -^2) — (C X yxy C X yyx Cy X xy Cy X yx) /2 
C{ E\) — {C X yxy ~\~ C X yyx ~\~ Cyxxy ~\~ Cyxyx) /2 

C{ E2) = (CxXXX Cxxyy Cyyxx ~\~ Cyyyy) /2 

C( ^4l) = (CxXXX ~\~ Cxxyy ~\~ Cyyxx ~\~ Cyyyy) /2, (3) 

where x, y correspond to e x , e y states. From this set of equations we find that the spacing between 
the singlets 1 A\ and 1 E 2 is equal to the spacing between the singlet 1 E\ and the ground state 
3 A 2 , i.e., CQ-Ax) - C( l E 2 ) = C( 1 £'i) - C( 3 A 2 ) = C XX yy + C yyxx = 2e, where the difference is 
the exchange energy. In addition, as 1 Ex and 1 E 2 belong to the same IR E, it can be shown 
that Ci}E 2 ) = C(}Ex) (see Appendix [bJ. Under this consideration, the ordering of the states is 
{ 3 A 2 , 1 E, 1 Ai} with relative energies {0,2e,4e}. It should be noted that, in this case, the most 
symmetric state has higher energy since the Coulomb interaction between two electrons is repulsive. 
This picture might be modified by the following effect. Since the Coulomb interaction transforms 
as the totally symmetric IR, the matrix elements between states with the same symmetry are 
non-zero. The states 1 E(e 2 ) and 1 E(ae) can couple via the Coulomb interaction, increasing the 
gap between them. A similar effect happens with the states 1 Ai(e 2 ) and 1 Ai(a 2 ). In Eq. ^ we 
did not take into account the effect of the other electrons present in the system. Nevertheless, our 
basic results here serve as a qualitative estimate for the energy of levels and provides useful insight 
into the structure of the NV center. The results of a very recent calculations based on many-body 
perturbation theory (MBPT) p3] supports our conclusion. 



IV. SPIN-ORBIT INTERACTION 



In the previous section the electronic spin did not directly enter into our considerations. For 
instance, the energy of the raj = 0, ±1 sublevels of the 3 A 2 ground state would have exactly the 
same energy. However, if the electronic spin is taken into account, one can infer from Table [I] 



s 



that in general the ms = and ms = ±1 projections transform as functions of different IRs. 
For example, in the ground state 3 A2, the ms = projection transforms as the IR Ax, while the 
ms = ±1 projections transform as the IR E. This implies that the projections do not share the 
same eigenenergies of the system. The spin-spin and spin-orbit interactions may result in splitting 
of these orbitally degenerate states. 

The spin-orbit interaction lifts the degeneracy of multiplets that have non-zero angular momen- 
tum, and is also responsible for transitions between terms with different spin states |21j . It is a 
relativistic effect due to the relative motion between electrons and nuclei. In the reference frame 
of the electron, the nuclear potential, <fi, produces a magnetic field equal to V(j) x v/c 2 . In SI units, 
this interaction is given by 

ff -=L4 (wxp) ©- (4) 

where V = ecf> is the nuclear potential energy, m e is the electron mass and p is the momentum. 
The presence of the crystal field breaks the rotational symmetry of this interaction. Since <j) is 
produced by the nuclear potential, it transforms as the totally symmetric representation Ax, and 
therefore W = (V x , V y , V z ) transforms as a vector, where Vi = dV/dxi. Since p also transforms as 
a vector, it is possible to identify the IRs to which the orbital operator components O = W x p = 
(VyPz — V z p y , V z p x — V x p z , V x p y — VyPx) belong. In C% v , the components of VV and p transform 
as (Ex, E%, A\) and therefore O transforms as the IRs (Ei, Ex, A2) = (E, The non-zero 

matrix elements of the orbital operators Oi in the basis {a, e x , e y } can be determined by checking 
if (ipi,Ok,(pf) D Ax and are shown in Table [TT| where A = (e y \O x \a) and B = (e x \O z \e y ) (for 
simplicity we denote by a the ai(2) orbital state ). In this case, the spin-orbit interaction can be 
written in terms of the angular momentum operators U and takes the following form: 

HsO = Ky (l X S x + lySy) + Kl Z S z , (5) 

where X X y(X z ) denotes the non-axial (axial) strength of the interaction. In a system with or 
spherical symmetry, A = B and the usual form (S • L) of the spin-orbit interaction is recovered. 
It is also useful to think about e± as p± orbitals and ai(2) as a p z orbital, where the angular 
momentum operators satisfy l±ax(2) oc e± [2^] . 

Once it is known how the spin-orbit interaction acts on the orbitals, e x , e y and a, it is possible 
to calculate the effect of this interaction on the 15 states given in Table [1} An important effect is 
the splitting in the excited state triplet between the states Ax,A2 and E x ,E y and between states 
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TABLE II: Matrix elements for orbital operators in the C^v point group. For the Td symmetry group or 
spherically symmetric potentials, A = B. 
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E x ,E y and E\,E2 [llj . The spin-orbit interaction can be written as, 

H SO = A„(|i4i}(Ai| + |A 2 )(A 2 | - - |^)(£? 2 |), (6) 

in the excited state triplet manifold {A±, A2, E x , E y , E±, E2} ■ Another effect, relevant when treating 
non-radiative transitions, is that the axial part of the spin-orbit interaction (A z ) links states with 
m s = spin projections among states of the same electronic configuration, while the non-axial 
part {\ x ,y) links states with non-zero spin projections with singlets among different electronic 
configurations. In Figure [T] we show the states linked by the axial and the non- axial parts of 
the spin-orbit interaction, for which non-radiative transitions might occur. In addition to the 
well known transition between A\(ae) — > 1 A\(e 2 ), we find that this interaction might also link 
E\-2 (ae) — > 1 E\^{e 2 ) and in particular E xy — > E x>y (ae). The latter transition may play an 
important role, as recent ab initio calculations have shown that the singlets 1 E x ^ y might lie very 
close in energy to the excited state triplet [H]. In our model, the non-axial part of the spin- 
orbit interaction, X xy (/+s_ + Ls+), does not mix the states of the excited state triplet with 
different spin projections because the raising and lower operators, /_ and /+, link states of different 
electronic configurations. In particular, this interaction cannot mix the states of the excited state 
triplet because the mixing is suppressed by the large energy gap that separates different electronic 
configurations. 

We have numerically evaluated the ratio between the axial part and transverse part of spin- 
orbit, X z /X xy = B/A = 0.75 using the functions e x and e y and oi(2) from ab initio calculations 
(see Appendix [E]) . This suggest that if the axial part of spin-orbit is 5.5 GHz [10], the non-axial 
part should be on the order of X xy = 7.3 GHz and only couples singlets with triplets states as 
shown in Figur^TJ We have also numerically confirmed the structure of Table [H] with three digits 
of precision in units of GHz (see Appendix ^1) . 
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V. SPIN-SPIN INTERACTION 



The spin-spin interaction between electrons is usually not present in systems with spherical 
symmetry, due to the traceless character of the magnetic dipole-dipole interaction. However, if the 
electron wavefunction is not spherically distributed, this interaction does not average out. Here 
we describe its effect on the excited state triplet of the NV center and we provide a numerical 
estimation of its strength. The spin-spin interaction can be written (in SI units) as, 

Mo 9 2 P 2 



(3(si • r)(s 2 • r) - si • s 2 ) , 



(7) 



Air r 3 

where Sj = | [a x , a y , a z ] are the spin operators of particle i and aj (j = x,y,z) are the Pauli 
matrices, (3 is the Bohr magneton, g is the Lande-factor for the electron and /io is the magnetic 
permeability of free space [27] . In order to analyze the effect of this interaction in the defect it is 
useful to write the spatial and spin parts separately in terms of the irreducible representations of 
the point group. Then, it is straightforward to express this interaction in terms of the eigenstates 
of the defect (see Appendix [C]) , 

H ss = &{\A 1 )(A l \ + \A 2 )(A 2 \ + \E l )(E l \ + \E 2 )(E 2 \) 
-2A(\E x )(E x \ + \E y )(E y \) 
+2&'{\A 2 ){A 2 \-\A 1 ){A 1 \) 

A"{\E l ){E y \ + \E y ){E l \-i\E 2 ){E x \+i\E x ){E 2 \), (8) 

where the gaps between the m s = ±1 and m s = projections and between A\ and A 2 states are 
given by 
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(9) 
(10) 



while the mixing term is given by 



A" 



3xz 
V^r 3 



X). 



(11) 



Figure [2] shows the effect of spin-orbit and spin-spin interactions on the excited state manifold. 
In particular, we find that the state A 2 has higher energy than the state A\ (2 A' > 0), contrary to 
previous estimations |11| [28] . In addition, we find that the spin-spin interaction A" mixes states 
with different spin-projections. This effect is the result of the lack of inversion symmetry of the NV 
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center and it is not present in systems with inversion symmetry such as free atoms or substitutional 
atoms in cubic lattices. This does not contradicts group theoretical estimates as the mixed states 
transform according to the same IR (e.g. the E\ and E y states both transform according to the 
IR Ei, see Table [J. 

We estimated these parameters using a simplified model consisting of the dangling bonds given 
in Figure [6] (in the Appendix) for the three carbons and the nitrogen atom around the vacancy. 
The dangling bonds are modeled by Gaussian orbitals that best fit to the wavefunction obtained by 
an ab initio DFT supercell calculation (see Appendix |E|), The distance between atoms is also taken 
from these simulations. To avoid numerical divergences when r = 0, we estimate Eq. ([9 11) in 



reciprocal space following Ref. [29]. The values for the zero field splitting (A es = 3A), gap between 
states Ai and A2 (4A') and mixing term between states E\^ and E x ^ y (A") are given in Figure [2jb. 
As ab initio calculations cannot accurately estimate the nitrogen population p^ = |/3| 2 in the single 
orbital state a\{2) (see Appendix [A] for a definition of parameter /3), we have plotted in Figure^, 
the values of the spin-spin interaction as a function of pn- In addition, the solid regions in the 
figure take into account variations of the relative distance among the three carbons, the nitrogen 
and the vacancy. The distance between the carbons and the vacancy is increased between and 
3%, meanwhile the distance between the nitrogen and the vacancy is decrease between and 4% 
relative to their excited state configuration (solid lines). This shows how the spin-spin interaction 
depend on the distance between the atoms. 

We emphasize that, contrary to the ground state of the NV center, the splitting between A\ and 
A2 in the excited state exists because the spin-orbit interaction mixes the spin and spatial parts. 
In fact, at high temperatures, where the spin-orbit interaction averages out [12], and if the spatial 



part is given by + |Y)(Y~|, it can be checked by looking at Eq. (C2) that only the zero 

field splitting A es survives from the electronic spin-spin interaction, as confirmed by experiments 
\V2\ I24j . In addition, the spin-orbit interaction in the excited state, Eq. ([6]), can be written as 
Hso = K\X){Y\ — \Y)(X\) <g> (|aa)(aa| — \/3(3)(f3f3\), which also vanishes if the spatial part is given 
by \X)(X\ + \Y)(Y\. 

VI. SELECTION RULES AND SPIN-PHOTON ENTANGLEMENT SCHEMES 

Group theory tells that transitions are dipole allowed if the matrix element contains the totally 
symmetric IR, {ipf\e\(pi) D A\. In the case of the NV center (Cs v ), the only non-zero matrix 
elements are (a\x ■ r\e x ) and (a\y ■ r\e y ), from which it is straighforward to calculate the selection 
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rules among the 15 eigenstates given in Table |T] for the unperturbed center. This is shown in 



Table III These matrix elements have been confirmed by our first-principles calculations of these 
matrix elements in the velocity representation as well as by other authors only for the triplet 
transition [30]. In addition to the well known triplet-triplet transition [31 j . transitions are allowed 
between singlets of different electronic configurations. We remark that the transition between 
singlet 1 A\(e 2 ) and singlet 1 E(e 2 ) is not strictly forbidden by group theory to first order, but since 
both states belong to the same electronic configuration, no dipole moment exists between them 
and the probability of radiative transition is extremely low. According to our results using wave 
functions from first-principles calculations (see Appendix [E]), the ratio between the dipole transition 
matrix elements associated with the singlet states to those of the triplet states is about 5xl0 -9 . 
The singlet-singlet transition might be allowed by phonons or mixing of the states with singlets 
of different electronic configurations. Recent experiments by Rogers et al. identified an emission 
from singlet to singlet [T2]. which we suggest is related to the 1 E(ae) — t- 1 A\(e 2 ) transition. The 
transition 1 A\{e 2 ) — y 1 E(e 2 ) might be possible for the reasons described above, but it is unlikely 
to be sizable. A recent MBPT calculation supports our conclusion [14J. A suitable experiment 
to unravel this issue would be to look at the presence of this emission under resonant excitation. 
In this case, if the state 1 E(ae) is above the excited state triplet, the state l E(ae) will be hardly 
populated and therefore no singlet-singlet transition should be observed. 

Once the selection rules are known for the defect, it is possible to realize interesting applications 
such as spin-photon entanglement generation |32j . In the case of the NV center, the system can 
be prepared in the ^(ae) state. Next, the electron can spontaneously decay to the ground state 
3 ^2- by emitting a photon with a + (right circular) polarization or to the state 3 ^2+ by emitting 
a (7_ polarized photon (see Figure [3]). As a result, the spin of the electron is entangled with the 
polarization (spin) of the photon. The implementation of this scheme is sensitive to strain, which 



will be analyzed in Section VII However, in Section VIII, we recognize that the application of an 
electric field can be used to overcome some of these issues and facilitate the next step of entangling 
between two NV centers. 

VII. THE EFFECT OF STRAIN 

Strain refers to the displacement An of the atomic positions when the crystal is stretched stretch 
Ax |26) . It is a dimensionless tensor expressing the fractional change under stretching, e^- = 
and it can be produced by stress (forces applied to the solid structure), electric field, or temperature 
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TABLE III: Selection rules for optical transitions between: the triplet excited state (ae) and the triplet 
ground state (e 2 ), the singlets (ae) and the singlets (e 2 ), and the singlet (a 2 ) and the singlets (ae). Linear 
polarizations are represented by x and y, while circular polarizations are represented by a± — x ± iy. As 
an example, a photon with <j + polarization is emitted when the electron decays from state A 2 (ae) to state 
3 ^ 2 -(e 2 ). 
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|33j . A systematic study of strain can be used to unravel the symmetry of defects and explore their 
properties |34| . Strain can shift the energy of the states as well as mix them. It can reduce the 
symmetry of the crystal field by displacing the atoms. However, not all nine components of strain 
change the defect in a noticeable way. The antisymmetric part of e« transforms as a generator of 
the rotational group and therefore only rotates the whole structure. The symmetry and energies of 
the unperturbed states do not change upon rotation. Only the symmetric part of strain, e = e + e T 
affect the structure of a defect [26]. As with any other element of the theory, strain can be expressed 
in terms of matrices that transform according to the IRs of the point group under consideration. 
These matrices can be found by projecting a general strain matrix on each IR, 

e r = ^ X :^e. (12) 

e 

In Appendix [D] we show in detail how to deteermine the effect of strain on the eigenstates of the 
defect. For simplicity, in the case of the NV center we only write the effect of strain in the manifold 

Hstrain = ^Al^l + ^Al^l + ^El^l + ^E2^2 + ^El^l + ^E2^2 (13) 

where 5 a Al = (e xx + e yy )/2, 5 b A1 = e zz , 5% x = (e xx -e yy )/2, b% 2 = (e xy + e yx )/2, 5 b E1 = (e xz + e zx )/2, 
S E2 = ( e yz + e zy )/2 and 

(1 
1 


(000 

1 

in the manifold {e x ,e y ,a}. The effect of strain on the orbitals a,e x ,e y is easy to see. Af will 
shift equally the energies of the states e x and e y , while A\ will shift the energy of states a. Note 
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that both describe axial stress: the former leaves the e 2 electronic configuration unaffected and 
the latter leaves the ea configuration unaffected. Either one produces relative shifts between both 
configurations, resulting in an inhomogeneous broadening of the optical transitions. However, they 
do not change the selection rules. Only the stress Af + A\, corresponding to either expansion or 
contraction, leaves all relative energies unaffected. E% splits the energy between e x and e y and Ey 
mixes the two states. Finally, E\ and Ey mixes the states e x and a± and e y and a\, respectively. 
In the case of the NV center, the effect of the matrices E x y can be neglected thanks to the large 
gap between orbitals a and e x ,y Therefore, in what follows we do not consider them further. 

Recent work has been done to analyze how strain affects the excited state structure of the 
NV center \10\ I12j . Here we derive the explicit form of strain affecting the different electronic 



configurations and look at how strain affects the selection rules described in Section VI 



The relevant strain matrices we will consider are E% and for which the Hamiltonian is, 



Hgtrain = S E1 (l e a;)( e a:| - \e y ){e y \) + 5 E2 (\e x )(e y \ + \e y ) (e x \) . 
This mostly affects the singlet and excited state configurations in the following form, 



(15) 
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(16) 



/ 

for the manifolds {A%, A 2 , E x , E y , E x , E 2 }, { 1 E 1 , 1 E 2 , 1 Ai} and ^Ex^Ey}, respectively. The 
ground state, due to its antisymmetric combination between e x and e y , is stable under the pertur- 
bation H s train- This can be checked by applying Eq. (15) to the ground state given in Table |l} 



The effect on the excited state triplet can be seen in Figure |4^, where the unperturbed states are 
mixed in such a way that, in the limit of high strain, the excited triplet structure splits into two 
triplets with spatial wavefunctions E x and E y . When strain overcomes the spin-orbit interaction 
(S%i > 5.5 GHz), the spin part decouples from the spatial part and the total angular momentum 
is no longer a good quantum number. Transitions from the excited state triplet to the ground 
state triplet are linearly polarized, where the polarization indicates the direction of strain in the 
xy plane. 

Figure Hfc shows how the polarization of the emitted photon from the state A 2 to the ground 
state 3 vl2-) varies from circular to linear as a function of strain. In the case of 5 E2 strain, the 
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effect is similar but now the mixing is different. As shown in Figure [4j A2 mixes with E\ and 
the photons become polarized along x — y. Note that, in the limit of low strain, in both cases the 
polarization remains right circularly polarized for the transition between the excited state ^(ae) 
to the ground state 3 A2- (e 2 ), while the polarization remains left circular for the transition between 
the excited state A2(ae) to the ground state 3 A2+(e 2 ). The fact that at lower strain the character 
of the polarization remains circular has been succesfully used in entanglement schemes [7j- The 
polarization properties of the states Ex 2 are similar to those of the states Ax, 2 but with the opposite 
polarization. 



VIII. STRAIN AND ELECTRIC FIELD 



The application of an electric field to a defect leads to two main effects. The first effect, the 
electronic effect, consists of the polarization of the electron cloud of the defect, and the second 
one, the ionic effect, consists of the relative motion of the ions. It has been shown that the two 
effects are indistinguishable, as they have the same symmetry properties [35]. The ionic effect is 
related to the well-known piezoelectric effect. When a crystal is under stress, a net polarization 
Pi = dijk&jk is induced inside the crystal, where d^ is the third-rank piezoelectric tensor and 
o-jk represents the magnitude and direction of the applied force. Conversely, the application of an 
electric field might induce strain given by ejk = d^Ei, where E{ are the components of the electric 
field [33 J . The tensor d^ transforms as the coordinates XiXjXj. and, therefore, group theory can 
be used to establish relations between its components for a given point group. In particular, the 
non-zero components should transform as the irreducible representation A\. By projecting dijk (or 
XiXjXk) onto the irreducible representation Ax, we can determine the non-zero free parameters of 
the tensor d and determine the effect of electric field on the eigenstates of the unperturbed defect 
(see Appendix [Dj . In the case of the NV center, the effect on the excited state triplet is given by 
following matrix, 
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in the basis {Ax, A%, E x , E y , E±, E2}, while the effect on the ground state triplet is 

H E = 2gbE z , (18) 

in the basis { 3 t42+, 3 A2o, 3 ^2-}- The parameters a, b and d are the components of the piezo electric 
tensor dijj. and g is the coupling between the strain tensor e and the NV center. Comparing Eq. 



(17) and (18), we note that the linear response of the excited state and ground state are in principle 
different. An electric field along the z (NV-axis) can be used to tune the optical transition without 
distorting the C% v symmetry of the defect, provided b 7^ d. In Figure [5^, we show the linear 
response of NV centers under an electric field parallel to the NV-axis. In this case, the linearity is 
not affected by the presence of strain. Our estimates for the ionic effect, based on the response of 
the lattice defect to electric field and the response of the orbital energies to strain (see Appendix 
[Pj ), indicate that the relative shift between the ground and excited state is about 4 GHz / MV/m. 
This could be very important in schemes to entangle two NV centers optically as the wavelength 
of the photons emitted from each NV center need to overlap [36] ■ In addition, an electric field with 
components E x ,y can be used to completely restore the C% v character of the defect. In Figure [SJa, 
we show the response of optical transitions under an electric field perpendicular to the NV axis. In 
this case, the response is linear if strain is absent and quadratic if strain is non-zero. Dashed lines 
show the response to an electric field when the defect experiences a 0.3 GHz strain along the [01-1] 
axis. Our estimations can be used to interpret the Stark shift observations by Tamarat et al. |37| . 



IX. CONCLUSIONS 



We have used group theory to identify, analyze and predict the properties of NV centers in 
diamond. This analysis can be extended to other deep defects in solids. A careful analysis of 
the properties of a defect using group theory is essential for predicting spin-photon entanglement 
generation and for controlling the properties of NV centers in the presence of perturbations such 
as undesired strain. We have shown that group theoretical approaches can be applied to determine 
the ordering of the singlets in the (e 2 ) electronic configuration and to understand the effect of 
spin-orbit, spin-spin and strain interactions. 
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Appendix A: Dangling bond representation and character table 

In this appendix we show in detail how to find the electronic representation for the case of the 
NV center. The NV center contains a vacancy that results in broken bonds in the system. In 
the tight binding picture, this means that three C atoms and one N-atom do not have enough 
immediate neighbor atoms to form a covalent bond for each of their valence electrons. These 
unpaired electrons are called 'dangling bonds'. In the case of the NV center, we consider a simple 
model consisting of four sp 3 dangling bonds, where three of them are centered on each of the three 
carbon atoms around the vacancy and the fourth dangling bond is associated with the nitrogen 
atom. The point group symmetry is C% v and its elements are the identity, rotations around the 
z (NV-axis) by ±27r/3 and three vertical reflection planes where each contains one of the carbons 
and the nitrogen. 

As discussed in Section |IIJ it is possible to construct the representation of the dangling bonds 
for the point group they belong to. Consider Figure [6] where the z axis is pointing out of the 
paper. The dangling bonds {01, 02, 03, a^} transform into one another under the operations of the 
C% v group. In this representation, each operation can be written as a 4 by 4 matrix, as shown in 
Figure |6| As representations depend on the particular choice of basis, it is customary to designate 
them using the trace of each matrix (characters). Note that the character for matrices belonging 
to the same class is the same, so in short the character representation for the dangling bonds is 
T a = {412}. This representation is clearly reducible, as it can be decomposed by the irreducible 



representation of the C% v group given in Table IV [38 



Application of Eq. ([!]) gives the following combination of ex's: {ac = (01 + 02 + 03) /3, e x = 
(201 — 02 — 03) /Vg, e y = (02 — 03) V% Qjv = ^n}, where ac and a^v transform as the totally 
symmetric irreducible representation A±, and e x and e y transform as functions of the IR E. Note 
that the e states transform as vectors in the plane perpendicular to the NV axis. 

Next, we model the electron-ion interaction to find out the ordering of these states. This 
interaction can be written in the basis of the dangling bonds 0, as, 

V = t>n|0Ar)(0jv| + ^2 l v i \a i ){a i \ + h n \ai)(a N \ + \ai){aj\h c (Al) 

i i>j 
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TABLE IV: Character and bases table for the double C 3v group. Examples of functions that transform 
under a particular representation are {z,x 2 + y 2 , z 2 } , which transform as the IR A\, the rotation operator 
R z as A 2l and the pair of functions {(x,y),(R x ,R y ),(xy,x 2 — y 2 ),(yz,xz)} as E. The spin projections 
{a(x),/3(i.)} transform as the IR Exj 2 (or -D1/2), while the functions aaa + i(3(3[3 and aaa — if3/3/3 transform 
as the IRs 1 E 3 / 2 and 2 E 3 / 2 , respectively. 
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where V{ < is the Coulomb interaction of orbital cij at site i, h c is the expectation value of the 
interaction between orbitals <7j and <7j + i at site i = {1, 2, 3}, v n = {ctn\ V\&n) and h n = {(Ti\V\a^f). 
This interaction, which transforms as the totally symmetric IR A%, not only sets the order of 
the orbitals but also mixes orbitals and ac- This is a consequence of the important concept 
that whenever a matrix element contains the totally symmetric representation, its expectation 
value might be different from zero [16]. Since both wave functions as well as the interaction 
between them transform as the totally symmetric representation At, the representation for the 
matrix element also transform as A%: Tq = T a (8) T aN Ti n t = A\ D A\. This interaction 
leads to the new basis[19] {a\{l) = aa c + (3a n , a\{2) = aa n + /3a c , e x = (2<n — o"2 — (73) /v6, 
e y = (o" 2 - ct 3 ) \/2}, with energies {^ai(i),oi(2) = \{v c +^h c +v n )±\/^, v c -h c , v c -h c }, respectively, 
where A = a/ (v c + 2h c — v n ) 2 + 12h"^, a 2 = 1—0 2 = 3/t^/A-E ai m. We see that the most symmetric 
state is lowest in energy, which is usually the case for attractive interactions. 

Appendix B: Ordering of singlet states 

Here we show that two states belonging to the same irreducible representation should have the 
same expectation value for their Coulomb interaction. We first note that the expectation value 
of an operator is a scalar and it should not depend on the particular coordinate system in use. 
In particular, this expectation value should be invariant under any operation of the Cs v group of 
the NV center. The Coulomb interaction is totally symmetric and therefore not affected by any 
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rotation, and the wavefunctions {e x ,e y } transform as the irreducible representation E. Therefore, 
we can get more information about these expectation values by projecting them on the totally 
symmetric irreducible representation A\, 

I h 

(ab, V, cd) = -J2 Xe (P R (a)P R (b), V, P R (c)P R (d)) . (Bl) 

II R=l 

We find as expected that 

CEx, V, 1 E 1 ) = -fEu V, X E X ) + ^E 2 , V, l E 2 ), (B2) 

which means that the states (^i, V, ^i) and ( 1 E 2 , V, 1 E 2 ) have the same energy, as required by 
symmetry. 



Appendix C: Spin-spin interaction 

In order to analyze the effect of spin-spin interactions (Eq. [7]) from the perspective of group 
theory, we first rewrite this interaction to identify spatial and spin terms that transform as IR 
objects in the point group, 

-3z 2 



h MQff 2 /3 2 
4tt 



-{si+s 2 - + si_s 2 + - 4:Si z s 2z 



3 x 2 - y 2 . 

+ | r 3 ( S l- fi 2- + S1+S2+) 

. 3 xy . 
+«^-5-(si-S 2 - - si + s 2+ ) 

3 xz 

+ n^( S l- S 22 + Sl z S 2 - + S 1 + S 2z + S lz S 2+ ) 

.3yz. . 

+' l n^{ S l- S 2z + Sl z S 2 - - Sl + S 2z - Sl z S 2+ ) 

where x, y and z are directional cosines and s± = s x ± is y . In the case of C$ v , for the unperturbed 
center, the expectation values of the 4th and 5th terms are nonzero in the spatial manifold of 
the excited state because the center lacks inversion symmetry. However, these terms 

might be neglected when considering other defects with inversion symmetry. We note now that 
the spatial part of the first term transforms as the totally symmetric representation A\, while the 
2nd and 3rd terms transform as the irreducible representation E. The reader can check which IR 
these combinations belong to by looking at the character table in the Appendix [A} Therefore, their 
expectation values can be written as 

^l 1 -^-) =A(\X)(X\ + \Y)(Y\) 
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(^#) = A'(|X)(X| - \Y)(Y\) (CI) 

£^ 2 (^^)-A'(|x)(y| + |y) ( x|), 
^^( m ^ m ) = A"(|r><Y|-|X)<X|), 

where |X) and |Y} are the two electron states given in Table|TJ Note that, for symmetry reasons, the 
second and third relations are characterized by the same parameter A', while the last two relations 
are characterized by the same parameter A". Similarly, it is possible to write the spin operators in 
the spin basis of the two holes, {|cra), \af3), \/3a), \/3/3)}. For example, S1+S2- = \a(3){/3a\. Using 



these relations and Eq. (CI ), the Hamiltonian in the fundamental bases of the excited state of the 
NV center is 

H ss = -A(\X)(X\ + \Y)(Y\) 

<g> {\aa)(aa\ + \PP){PP\ - 2\a/3 + pa)(ap + /3a\) 
-A'(\X)(X\ - \Y)(Y\) ® (\aa)(f3f3\ + \f3f3)(aa\) 
-iA'(\X){Y\ + \Y){X\) (\PP)(aa\ - \aa){P0\) 
+A"(\Y)(Y\-\X)(X\) 

® (\af3 + (3a)(aa - /3/3\ + \aa - Pf3){aj3 + (3a\) 
+iA"(\Y)(Y\-\X)(X\) 

8) (|a/3 + /3a) (aa + pp\ - \aa + PP)(aP + Pa\) . (C2) 

Finally, we can write H ss in terms of the eigenstates of the unperturbed defect (see Table [I]) . This 
leads to Eq. @. 

Appendix D: Strain and electric field 

The effect of strain on the electronic structure of the defect can be obtained from the effect 
of the electron- nuclei Coulomb interaction on the eigenstates of the defect. In our example, the 



Coulomb interaction is given by Eq. ( Al ). However, when the positions of the atoms are such that 
the symmetry of the defect is reduced, we should allow for different expectation values of the matrix 
elements: hij = {&i\V\(Tj) and hi n = {(Ti\V\(tn). We have assumed that the self interactions, v c and 
v n , do not change as the electrons follow the position of the ion according to the Born Oppenheimer 
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approximation. To relate the matrix elements to the ionic displacements, we can assume as a first 
approximation that the electron orbitals are spherical functions, and therefore the matrix elements 
can be parametrized by the distance between ions, hij qj) = hij (\q{ — qj\), so that we can write 



hij (|%'|) ~ hj 



ij 



+ 



1 dh 



{li ~ Qj) 



(5qi - 5qj) + .. 



(Dl) 



\Qij\ dqij 

The change in the matrix elements is linear in the atomic displacements. In turn, the atomic 
displacements are related to the strain tensor by 5qi = eqi, and therefore the change in the matrix 
element is given by 



Shij (\qij\) » 1 ^ {qi - qjf e (g 4 - qj) 



(D2) 



Under these considerations, it is straightforward to calculate the effect of strain on the eigenstates 
of the defect. For simplicity, we write here only the effect of strain on the degenerate orbitals, e x 
and e y , 

( 

5V = -g 



(D3) 



V e xy e yy 



where g = ^§|p- and q is the nearest neighbor distance between atoms. Using the electron 
wavefunction obtained from ab initio calculations (see Appendix [E]) we estimate that g ~ 2 PHz 
(P = peta = 10 15 ). 

The effect of electric field on the eigenstates of the defect can be analyzed by the inverse 



piezoelectric effect as described in Section VIII In this appendix we show how group theory can 
tell us the nature of the piezoelectric tensor. By projecting dijk (or XiXjXk) onto the irreducible 
representation A\, we can build the following relations, 



a — dm — —di2\ — 

b = diis = C?223 



-dl22 d — C?333 

C = dm = G?232 



(D4) 
(D5) 



and the d tensor can be written in the following short notation (contracted matrix form) [33] 



di 



ijk 



a —a 
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(D6) 



For a given electric field, we have a strain tensor of the form 



aE x + bE z 
—aE,, 



— a,Ey cE x 
-aE x + bE z cE y 



\ cE x 



CE; 
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dE z 



(D7) 
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To evaluate the magnitude of the piezo-electric response, we have used first-principles calculations 
as described in Appendix [Ej The values for the components of the piezo-electric tensor due to ionic 
effect are a ra b ra c ra 0.3 ^(MV/m)- 1 and d ra 3 ^(MV/m) -1 . 

Appendix E: Information about the first principles methods applied in our study 

To determine the values of the constants a,b,c and d introduced in Appendix |Dj we applied 
density functional theory (DFT) [39J calculations within a generalized gradient approximation PBE 
(Perdew-Burke-Ernzerhof) |40j . In the study of spin-orbit and spin-spin interactions we used a 512- 
atom supercell to model the negatively charged nitrogen-vacancy defect in diamond. Particularly, 
we utilized the VASP code |411l42j to determine the geometry of the defect which uses the projector 
augmented wave method [43} |4"4"| to eliminate the core electrons, while a plane wave basis set is 
employed for expanding the valence wavefunctions. We applied the standard VASP projectors for 
the carbon and nitrogen atoms with a plane wave cut-off of 420 eV. The geometry optimization was 
stopped when the magnitude of the forces on the atoms was lower than 0.01 eV/A. We calculated 
the geometry of both ground and excited states. We applied the constrained DFT method to 
calculate the charge density of the excited state, that is, by promoting one electron from the a\{2) 
orbital to the e x ,e y orbitals as explained in Refs. [251 US]- This procedure is a relatively good 
approximation as confirmed by a recent many-body perturbation theory study [H]. The obtained 
geometries from VASP calculations were used as starting points in the calculations of spin-orbit 
and spin-spin interactions. 

The spin-orbit energy was calculated by following Eq. Q in our manuscript. Since the spin-orbit 
interaction is short-range, we applied all-electron methods beyond the frozen-core approximation. 
We utilized the CRYSTAL code [46J for this calculation using the PBE functional within DFT. We 
took the geometry as obtained from the VASP calculation. We applied 6-31*G Gaussian basis set 
for both the carbon and nitrogen atoms. The calculated properties (like the position of the defect 
levels in the gap) agreed well with those from plane- wave calculations. We obtained the all-electron 
single particle states and the corresponding Kohn-Sham potentials on a grid and calculated the 
spin-orbit energy numerically. 

Finally, we also studied the piezo-electric effect. In this case an external electric field was 
applied along the NV-axis and perpendicular to it. For this investigation only a finite size model 
can be used, thus we modeled the negatively charged nitrogen- vacancy defect in a molecular cluster 
consisting of 70 carbon atoms and one nitrogen atom. The defect was placed in the middle of the 
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cluster. The surface dangling bonds of the cluster were terminated by hydrogen atoms. In our 
previous studies we showed [19] that the defect wave functions are strongly localized around the 
core of the defect, thus our cluster model can describe reasonably well the situation occuring in 
the bulk environment. For this investigation we again applied DFT with the PBE functional as 
implemented in the SIESTA code |47j . We used the standard double-^ polarized basis set and 
Troullier-Martins norm-conserving pseudopotentials [15J. This method gives identical results with 
those obtained from plane wave calculations regarding the geometry and the wave functions in 
supercell models [19] . We fully optimized the defective nanodiamond with and without the applied 
electric field. In this case we applied a very strict limit to the maximum magnitude of forces on the 
atoms, 0.005 eV/A. We applied 6 different values of the external electric field along the NV-axis and 
in perpendicular directions to it, where we could clearly detect the slope of the curvature of atomic 
displacements versus the applied electric field. The resulting values for the atom displacements in 
the presence of 1 MV/m electric field are on the order of a few 0.1 yuA. 
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FIG. 1: Energy diagram of the unperturbed nitrogen- vacancy center in diamond. Note that 
each electronic configuration can contain triplets (left column) as well as singlets (right column) which have 
been drawn in separated columns for clarity. Red arrows indicate allowed optical transitions via electric 
dipole moment interactions. The circular arrows between the states E\p. and E XjV represent the mixing 
due to spin-spin interaction (see Figure [2]). Dashed lines indicate possible non-radiate processes assisted by 
spin-orbit interaction. In the ground state (e 2 configuration), the distance between singlets and triplets is 
equal to the exchange energy of Coulomb interaction (2e). The horizontal dashed blue line represents the 
orbital energy of the ground state (without including spin-spin interaction). 
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FIG. 2: Splitting due to spin-orbit and spin-spin interaction in triplet ae. (a) The axial part of 
the spin-orbit interaction splits the states {Ai, A 2 }, {E x , E y } and {E\, E2} by A z . The spin-spin interaction 
splits states with different spin projections and also splits the A\ and A 2 states. Our theory predicts the 
A2 state at higher energy than the Ai state and that the states (£1,2) and (E XtV ) are mixed. As the state 
A\ has an additional non-radiative decay channel, it is possible to confirm this finding by measuring the 
lifetime of the state. Note that the splitting between A\ and A 2 is a direct consequence of spin-orbit mixing 
the spatial and spin part of the wavefunction. (b) Values for the zero field splitting (3A), gap between 
the states A\ and A 2 (4A') and mixing term (A") due to spin-spin interaction in the excited state as a 
function of the nitrogen population, p^, in the state oi(2). The shadowed areas indicate the possible values 
for these parameters when the distance between the vacancy and the three carbons is increased between 
and 3%, and the distance between the vacancy and the nitrogen is decreased between and 4% of their 
excited state configuration. The solid lines correspond to the maximum (minimum) distance between the 
carbons (nitrogen) and the vacancy. 
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FIG. 3: Spin-photon entanglement generation. When the NV center is prepared in the excited state 
A2( 3 E), the electron can decay to the ground state 3 ' Ai m s = 1 [m s = —1) by emitting a right (left) 
circularly polarized photon. 
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FIG. 4: Excited state structure as a function of strain, (a) Eigenvalues of the excited state triplet as 
a function of 8% x strain, (b) Mixture of the eigenstate with higher energy (corresponding to A2 in the limit 
of low strain) and (c) the polarization of dipolar radiation under transitions from this state to the 3 ^2+ 
state of the ground state. Note that in both cases the circular polarization character of radiation remains. 
On the other hand, the linear polarization rotates 90° for strain along 5 b E2 with respect to that of strain 
along 8% 1 . 
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FIG. 5: Piezo-electric response of optical transitions, (a) response to electric field E z along the NV- 
axis ([111] orientation or equivalents). The defect only shows linear Stark Shift independent on the initial 
strain, (b) electric field E x applied perpendicular to the NV-axis in the absence of strain (solid lines). The 
optical transitions 3 A 2 (m s = 0) — > E x {m s — 0) and 3 A 2 (m s = 0) — > E y (m s — 0) are split linearly and 
evenly. In the presence of strain along the y direction (dashed lines), the response is quadratic due to the 
splitting between E x and E y states in the excited state. Our numerical results are in fair agreement with 
experimental results [37]. 
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FIG. 6: Schematic of the NV defect and dangling bond representation. (Top) Schematics of the 
dangling bond orbitals used to represent the NV defect. The symmetry axis or NV axis is pointing out 
of the plane of the page. The dashed lines represent the three vertical reflections planes of the C 3v group. 
(Bottom) Matrix representation of the dangling bonds. 



